.  AD;-753 


3  58' 


SUPERSONIC  LAMINAR  FLOW  REATTACHMENT 
AND  REDEVELOPMENT  &EHIND  A'  TWO-DIMEN¬ 
SIONAL  REARWARD  FACING  STEP 

Donald  J.  Spring 

Army  Missile  Command 
Redstone  Arsenal,  Alabama 


31  October  1972 


i 


! 


I 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  of  commerce  , 

5285  Port  Royal  Road,  Springfield  Va.  22151 


i 


TECHNICAL  REPORT 
RD-73-16 


SUPERSONIC  LAMINAR  FLOW  REATTACHMENT 
AND  REDEVELOPMENT  BEHIND  A  TWO-DIMENSIONAL 
REARWARD  FACING  STEP 


Donald  J.  Spring 


October  1972 


f^~n,‘~rr!-v 

i-J  \~r  >-*•  % 

T^EJSrlKUE  [ 

^  DEC  26  1572 

I  h  ^iTQTlJTE  _r. 


\A^ 


Approved  for  public  release;  distribution  unlimited- 


mqssoile  €©&»» 

Fted&ton&  ^rserteil ,  A/abama 


NATIONAL  TECHNICAL 
INFORMATION  service 

O  S  C‘  r o~'.'Ttrrc* 

Sw'  -v’  <*  "  *  A  22151 


DEC  65  PREVINS  EDIT 5 ON  IS  OBSOLETE 


STS 

see 

t  inifiFi'Mliat  ..„ 


L-!s  “”«•»  □ 


i  . — 

f  «/R*^liA3iS,|ff  «?» 

j — — _ "AlL'tSrtTv  ?.n— 

I ./) !  T^; 


i  i 

1/  i 


DISPOSITION  INSTRUCTIONS 


DESTROY  THIS  REPORT  WHEN  IT  IS  NO  LONGER  NEEDED.  DO  NOT 
RETURN  IT  TO  THE  ORIGINATOR. 


DISCLAIMER 


THE  FINDINGS  >N  THIS  REPORT  ARE  NOT  TO  BE  CONSTRUED  AS  AN 
OFFICIAL  DEPARTMENT  OF  THE  ARMY  POSITION  UNLESS  SO  DESIG¬ 
NATED  BY  OTHER  AUTHORIZED  DOCUMENTS. 


TRADE  NAMES 


USE  OF  TRAOE  NAMES  OR  MANUFACTURERS  IN  THIS  REPORT  DOES 
NOT  CONSTITUTE  AN  OFFICIAL  INDORSEMENT  OR  APPROVAL  OF 
THE  USE  OF  SUCH  COMMERCIAL  HARDWARE  OR  SOFTWARE. 


I 


!  lb  1 1,1  i  1 1 1 1 1 1  d  fci  a  liiklilifiSBffl  0ttS33BiMAttiaili 


UNCLASSIFIED _  47 

DOCUMENT  CONTROL  DATA  •  R  &  D 

(Security  ct— Illation  of  t* tto,  body  of  rtiffcl  snd  indexing  onnotmtlon  mutt  bo  ontorod  when  tf>#  ov#r#7/  r«pwf 

II.  omCINATINQ  ACTIVITY  (Cotponi  Touthoe)  zT.  REPORT  SECURITY  CLASSIFICATION 

Aerobal lis tics  Directorate  UNCLASSIFIED 


US  Army  Missile  Res,  Dev  &  Eng  Lab 

US  Army  Missile  Command  1 

Redstone  Arsenal,  Alabama  35809 _ 

S.  REPORT  TITLE 

SUPERSONIC  LAMINAR  FLOW  REATTACHMENT  AND  REDEVELOPMENT 
BEHIND  A  TWO-DIMENSIONAL  REARWARD  FACING  STEP 

4.  descriptive  NOTES  (Typo  of  report  end  fneiuolvo  dmf») 

Technical  Report 

»•  «UTHOR(»l  (Flr.l  nun,  333!  Initial,  TSuTUSSSJ  ______________ 


Donald  J,  Spring 


■«.  REPORT  OATE 


1 7 A.  TOTAL  NO.  OF  PACES  1 76.  NO.  OF  REFS 


31  October  1972 

S*.  CONTRACT  OR  ORANT  NO. 


j^5~I 


M.  ORIOINATOR**  REPORT  NUWfBER(S) 


A.  PROJECT  NO. 


10.  DISTRIBUTION  STATEMENT 


RD-72-16 


•6.  OTHER  REPORT  NO(S>  (Any  otht  numboro  that  muty  bo  Ml/flFf 
Ifci#  roport) 


Approved  for  public  release;  distribution  unlimited. 


(  If.  SUPPLEMENTARY  NOTES 


1 12.  SPONSORING  MILITARY  ACTIVITY 


Same  as  No.  1. 


A  theoretical  study  of  the  flow  recompression,  reattachment,  and  redevelopment 
associated  with  a  supersonic  laminar  ieoenergetic  flow  past  a  bark  step  wau  made. 
Based  on  an  integral  approach,  the  governing  system  of  ordinary  differential 
equations  for  such  recompression  and  redevelopment  processes  was  derived  by  intro¬ 
ducing  simplified  flow  models.  The  preceding  constant  pressure  jet  mixing  region 
was  analyzed  so  as  to  provide  the  necessary  initial  conditions  lor  the  recompression 
process.  The  effect  of  finite  initial  boundary  layer  on  this  mixing  process  was 
-aleo  taken  into  account.  The  system  of  equations  for  recompression  was  integrated 
'with  the  base  pressure  and  the  starting  location  of  recompress { in  as  parameters. 

These  parameters  were  then  determined  through  the  conditions  to  be  satisfied  at  the 
point  of  reattachment. v  The  calculated  base  pressure  ratios  agreed  well  with  the 
experimental  data.  Some  analysis  was  continued  for  the  flow  i  development  downstream 
of  reattachment.  With  slight  modification  of  the  results  at  i.  e  point  cf  reattach¬ 
ment,  the  pressure  distribution  within  the  redeveloping  flow  region  can  be  predicted 
fairly  well  for  the  flow  case  of  Mj.  =  2.0.  Additional  possible  improvements  to  the 
analysis  was  suggested  and  discussed. 


-1473 


REPLACES  DO  FORM  Nil.  1  JAM  M,  WHICH  IS 
OBSOLETE  FOR  ARMY  USE. 


UNCLASSIFIED 

tecatity  QuSKoSm 


CONTENTS 


1.  INTRODUCTION  . 

2.  DEVELOPMENT  OF  ANALYTICAL  EQUATIONS  . 

2.1  General  Considerations  . 

2.2.1  Upper  Viscous  Layer  ...  . 

2.2.2  Wake  Flow  Region  . 

2.2.3  Laminar  Mixing  Analysis  . 

2.3  Analysis  for  Redevelopment  Region  . 

2.3.1  The  Flow  Condition  at  the  Initial  Section 

of  Redevelopment  .  .  . 

2.3.2  Downstream  Redevelopment  Flow  ..... 

3 .  CALCULATION  METHOD  AND  RESULTS  . 

3.1  Calculation  Procedures  . 

3.2  Results  of  Calculations  .  ,  . 

4.  DISCUSSION  AND  CONCLUSIONS  . 

REFERENCES  . 


iii 


Preceding  page  blank 


NOMENCLATURE 


(pu/peue),  Chapman-Rubesin  parameter 

V/V  ,  Crocco  number 
max 

Normalized  shear  function 

Step  height 

Height  for  back  flow 

Length  of  the  constant  pressure  region  along  the  wake 
boundary  measured  from  the  separation  corner 

H/sin  0  ,  length  of  the  wake  boundary  measured  from 
00 

the  separation  corner  to  the  point  of  reattachment 

Mach  number 

Pressure 

Specific  gas  constant 

U  /v  ,  Reynolds  number 
00  00 

Location  of  the  reattachment  point  measured  from  the 
base  of  the  step 

Streamwise  distance  function  defined  by  Eq.  (2.29) 

Temperature 

x-velocity  component 

Magnitude  of  velocity 

y-velocity  component 

Coordinate  in  main  flow  direction 

Coordinate  normal  to  x 

u  /u^;  reference  perturbation  velocity 


Streamline  angle 


Ratio  or  specific  heat 
Thickness  of  viscous  layer 
y/S  or  y/&b 
o^y/x,  mixing  ordinate 
Angle  or  momentum  thickness 
Viscosity 

Kinematic  viscosity 
Density 

Laminar  spread  rate  parameter 
Shear  stress 

u/u^,  dimensionless  velocity 
Prandtl-Meyer  function 


Viscous  layer  above  the  dividing  streamline 
Viscous  layer  below  the  dividing  streamline 
Dividing  streamline 
External  inviscid  stream 

Characteristic  length  for  Reynolds  number  measured 
from  leading  edge  of  plate  to  separation  point 

Wall 

Stagnation  state 

Approaching  flow  state,  or  initial  control  volume 
station 

Flow  state  after  the  Prandtl-Meyer  expansion,  or 
final  control  volume  station 

Station  at  the  beginning  of  recompression 
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1 .  INTRODUCTION 


One  of  the  problems  facing  designers  of  flight  vehicles  with 
relatively  blunt  bases  has  been  the  accurate  prediction  of  the  base 
pressure  acting  upon  the  body  base.  For  certain  configurations  and 
flight  conditions  the  resulting  base  drag  makes  a  significant  contri¬ 
bution  (as  much  as  one-third)  to  the  overall  drag  of  the  vehicle. 
Recognition  of  this  problem  has  led  to  a  considerable  amount  of 
research  over  the  last  quarter  of  a  century  in  an  attempt  to  reach  a 
better  understanding  of  the  flow  field. 

It  was  recognized  early  that,  for  a  flow  past  a  blunt  base  or 
step  (see  Fig.  1),  the  flow  is  unable  to  follow  the  body  contour,  thereby 
resulting  in  a  separation  of  the  flow.  For  supersonic  approaching 
flows,  it  will  expand  from  the  initial  pressure  to  that  of  the  base 
pressure  at  the  corner.  At  the  same  time,  the  boundary  layer  will 
also  follow  this  expansion  process.  As  a  result  of  the  separation,  the 
free  shear  layer  is  energized  by  a  mixing  process  and  prepares  itself 
for  the  subsequent  reattachment.  The  interaction  between  the  inviscid 
and  viscid  stream  controls  and  determines  the  overall  flow  field.  For 
practical  purposes,  the  pressure  in  the  initial  part  of  the  mixing 
region  is  constant ;  however,  close  to  the  recompression  region,  this 
simplification  is  not  justified  since  appreciable  pr-.  ".sure  rises  have 
been  observed  along  the  flow.  In  addition,  considerable  variation  of 
static  pressure  exists  across  the  flow  as  a  result  of  the  finite 
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curvature  of  the  streamlines.  Little  effort  has  been  expended  in 
understanding  the  flow  field  in  the  region  of  recompression  or 
reattachment. 

Early  researchers  [1,  2,  3]  developed  empirical  methods  to  pre¬ 
dict  the  base  pressure  which  was  found  as  only  a  function  of  free- 
stream  Mach  number  and  it  was  not  until  the  early  work  of  Chapman  [4] 
that  the  importance  of  Reynolds  number  effects  was  recognized. 

Chapman  proceeded  to  formulate  a  semiempirical  method  and  showed  the 
dependency  of  base  pressure  upon  Mach  number  and  characteristic 
Reynolds  number.  Further  extensions  of  these  ideas  were  proposed  by 
Love  [5].  These  semiempirical  methods,  while  giving  quick  estimations 
of  base  pressures,  were  unable  to  describe  the  mechanics  of  the  flow 
field. 

The  first  attempt  to  qualitatively  explain  the  complex  flow  field 
appears  to  have  been  accomplished  by  Crocco  and  Lees  [6] .  They  used 
an  integral  approach  where  the  momentum  transfer  across  the  mixing 
region  was  estimated  by  a  proper  average  constant  rate.  The  wake 
region  is  thus  energized  until  it  can  support  the  recompression  pro¬ 
cess  as  the  flow  returns  to  the  direction  of  the  original  flow. 

Chapman  [7]  suggested  the  concept  which  presently  is  the  basis  for  most 
of  the  component  approaches  to  the  base  pressure  problem.  He  visual¬ 
ized  a  region  at  the  base  of  the  body  commonly  referred  to  as  the 
"dead-air"  region  which  was  responsible  for  maintaining  equilibrium  in 
the  flow  and,  thus,  the  concept  of  a  dividing  streamline  was  born. 

Korst  [8j  adapted  the  dividing  streamline  concept  in  a  base  pressure 
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theory.  This  approach  is  conanonly  referred  to  as  the  Chapman-Korst 
model.  The  equilibrium  condition  states  that  the  total  pressure  along 
the  dividing  streamline  must  possess  only  enough  mechanical  energy  to 
support  the  recompression  process  up  to  the  rear  stagnation  point.  The 
fluid  beneath  the  dividing  streamline  is  returned  and  conserved  within 
the  wake  while  that  above  the  dividing  streamline  continues  downstream. 
Korst  used  an  eddy -viscosity  concept  to  solve  several  cases  of  turbu¬ 
lent  mixing  while  Chapman  [9]  used  a  similarity  solution  for  laminar 
mixing.  Both  methods  assumed  that  the  pressure  at  reattachment  could 
be  determined  from  i.sentropic  relationships. 

Following  these  efforts  Lees  et  al  [10,  11]  extended  the  Crocco- 
Lees  method  by  developing  a  technique  that  eliminated  the  experimen¬ 
tally  determined  mixing  coefficient  by  introducing  the  first  moment  of 
momentum  equation.  They  showed  that  a  throat  or  critical  point  existed 
in  the  flow  downstream  of  the  reattachment  point.  The  base  pressure 
was  determined  by  the  requirement  that  the  flow  pass  smoothly  through 
this  throat.  Alber  and  Lees  [12]  later  extended  the  integral  viscid- 
inviscid  theory  to  turbulent  flows.  A  turbulent  eddy-viscosity  model 
was  formulated  from  that  of  an  incompressible  flow  by  adopting  a  cer¬ 
tain  reference  density.  However,  in  a.  1  of  the  preceding  case,  no 
detailed  studies  of  the  recompression  process  were  completed.  Most 
of  the  other  methods  used  empirical  corrections  based  upon  experimental 
data  such  as  that  proposed  by  Nash  [13]  to  accurately  predict  the 
reattachraent  pressure . 

Later  efforts  directed  to  solving  the  wake  problem  have  included 
the  method  of  integral  relations  as  applied  by  Crawford  [14]  and 
numerical  solutions  to  the  Navier-Stokes  equations  [15] .  In  addition 
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Co  these  methods,  Weiss  and  Weinbaum  [16,  17,  18,  19]  obtained 
detailed  solutions  of  the  near  wake  by  dividing  the  flow  field  ’into 
three  regions.  Each  region  is  solved  using  a  different  technique, 
hence,  their  methods  are  commonly  referred  to  as  multimethod  base 'flow 
theories.  The  rotational  outer  flow  resulting  from  the  expansion  of 
the  boundary  layer  is  solved  using  the  method  of  characteristics, 
the  viscous  layer  above  the  dividing  streamline  is  treated  using  a 
modified  Oseen  solution  of  the  boundary  layer  equations,  and  the  recir¬ 
culation  region  is  solved  by  a  finite-difference  representation  of 

'  the  full  Navier-Stokes  equations.  Coupling  is  assured  by  requiring 

continuity  of  flow  properties  on  the  dividing  streamline.  The  disad¬ 
vantage  to  this  last  group  of  meth'  _■  is  that  even  with  the  advent  of 
high  speed  computers  a  considerable  amount  of  computational  time  is 

1 

I 

|  required  to  study  each  case. 

i  Chow  [20]  in  1570  developed  a  technique  to  study  the  turbulent 

i 

!  reattachment  which  yields  satisfactory  results  for  the  recompression 

pressure  distribution  up  to  the  reattachraent  point  without  the  exten¬ 
sive  use  of  empirical  information.  A  constant  pressure  turbulent  jet 
mixing  precedes  the  recompression  region  where  due  to  the  turning 

| 

1  process  the  normal  pressure  gradient  can  no  longer  be  neglected.  By 

\  developing  a  system  of  equations  which  also  must  be  satisfied  at  the 

reattachment  point, a  solution  can  be  obtained. 

The  present  effort  is  an  examination  of  the  reattachment  process 
within  tho  two-dimensional  laminar  flow  regime.  By  following,  in 
general,  the  approach  as  suggested  by  Chow  for  turbulent  flow,  it  is 
felt  that  laminar  flow  cases  can  be  studied  in  more  detail  since  many 
factors  which  influence  the  recompression  processes  and  had  to  be 
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ignored  in  turbulent  flow  studies,  can  be  investigated  and  explored. 

In  addition,  the  study  of  flow  redevelopment  and  rehabilitation  after 
reattachment  which  can  not  be  adequately  handled  in  turbulent  flows 
would  fully  illustrate  the  behavior  of  relaxation  toward  the  equilib¬ 
rium  flow  far  downstream.  Isoenergetic  flow  shall  be  assumed  through¬ 
out  this  analysis  and,  thus,  the  energy  consideration  is  conveniently 
eliminated. 
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2-  DEVELOPMENT  OF  ANALYTICAL  EQUATIONS 

•  i 

\ 

2.1  General  Considerations  : 

.  I 

It  is  expected,  that  the  boundary  layer  concepts  still  give  a 
valid  description  of  the  flow  for  ;he  present1  problem.  However,  it  is 
also  recognized  that  within  the  region  of  recompression  and  redevelop¬ 
ment,  the  curvature  of  the  streamlines  is  not  negligible;  i.e.,  the 

simplification  of  Sp/5y  =  0  across  the  viscous  region  can  no^longer 

? 

be  adopted.  Integral  analysis  was  employed  to  solve  the  governing 
system  of  equations  derived  for  the  viscous  layers  incorporating  their 

i 

interaction  with  the  outer  free  stream.  The  requirement  that  the  s:ame 
flow  properties, along  the  common  boundary  joining  the  upper  and  lower 
viscous  layers, determines  a  proper  solution .  .  > 

i 

2.2  Formulation  of  the  Jet  Reattachment  Process 

The  recompression  region  shown  in  Fig.  2  is  divided  into  two 
parts  along  the  dividing  streamline.  Fluid  contaihed  within  the  region 
above  the  dividing  streamline  will  proceed  downstream  while  that  below 
the  dividing'  streamline  will  be  turned  back  into  the  wake  region.  The 
upper  shear  layer  interacts  with  the  external  flow  which  guides  and 
receives  the  influence  of  the,  shear  layer  by  following  an  inverse 
Prandtl-Meyer  relationship.  One  can  describe  the  interaction  from  the 
fact  that  the  transverse  velocity  component  at  the  edge  of  the  shear 
layer  induces  an  increase  in  pressure  in  the  outer  free  stream  which 
in  turn  exerts  an  influence  on  the  flow  properties  within  the  layer  ‘ 
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including  this  transverse  velocity  component  at  the  edge  of  the  shear 
layer  itself. 

2.2.1  Upper  Viscous  Layer 

Conservation  of  mass  across  the  upper  shear  layer  (see  Fig.  3) 
requires 


"3  (ou)  +  ^  (ov)  =  0 
oa  uy 


(2.1) 


Choosing  an  x-coordinate  to  be  coincident  with  the  dividing  streamline, 
one  can  integrate  Eq.  (2.1)  across  the  shear  layer  to  obtain 


o  v 
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S  <°u)  dy 
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Upon  employing  Leibnitz's  Rule,  the  abo^e  equation  can  be  rewritten  as 


p  v 
Me  e 


rl 

=  -  / 
dx  J 


d& 


pu  dy  -  p  u 


e  e  dx 


I  . 


m  e 


After  nondimens ionalizing  with  the  free-stream  properties  and  after 
some  manipulations,  one  obtains  the  streamline  angle  at  the  edge  of 
the  upper  shear  layer  as  given  by 


-  o  e  d 

tan  3  =  —  =  — 

e  u  dx 
e 


,/h. 


dt 

e _ d_ 

1  dx 


k  -  0 


(7-1) 
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at  t,  =  0 


at  t;  =  1 


P  =  Pd 


p  =  pe; 


ve  / 


Eq.  (2.6)  can  now  be  written  as 


KT\\  f  r'?(1'  *>  d« +  Mk6.  / 


•  4(^)4 -3 (‘--3 


For  this  upper  shear  layer,  a  velocity  profile  of  third-order 
polynomials  is  assumed,  i.e., 

*p  *  ^ +  If  Id  5  +  3(x  -  «Pd)  - 2  if  |d  +  if  d  -  2(x  -  ^d)  >  ( 


which  satisfies  the  boundary  conditions 


'p =  ^  5  if =  if  |d  at  s  = 0  » 


(** "  i J  if " 0  at  s  = 1 

It  is  recognized  that  at  the  rear  stagnation  point  (<p^  =0),  th 
velocity  and  the  stress  should  vanish  so  that  a  simple  correlation 
between  the  slope  parameter  and  qjj  would  be  that  they  are 

linearly  proportional  to  each  other,  with  the  constant  of  proportion 
ality  being  evaluated  from  initial  conditions  prevailing  at  the  end 
of  the  mixing  region. 


11 


The  expression  for  the  pressure  difference  across  the  upper  layer 
can  be  obtained  from  vertical  momentum  considerations.  The  basic 
equation  is  given  by 


JL 

3x 


(puv)  +  ^ 


(2.9) 


Integrating  across  the  upper  layer  one  obtains 


_d_ 

dx 


A 

2,  /  p  2 

u  o  I  —  © 

e  aJ  Pe 


tan  p  d£ 


^a  2  2 

“puv  -z - h  p  u  tan  p 

Me  e  e  dx  Ke  e  e 


5 

'  (pd  '  pe)  +  S  f  pvd“  1 


(2.10) 


with  tan  p  =  tan  pe(2{;  -  £  2)  . 

The  shear  stress  in  Eq.  (2.9)  has  been  evaluated  from 


T 


xy 


where  u  is  the  viscosity  which  is  assumed  to  be  linearly  proportional 
to  temperature.  Upon  rearranging  and  normalizing,  Eq.  (2.10)  yields 
the  final  expression  for  the  pressure  difference  across  the  upper 
layer,  i.e., 


(2.11) 
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2.2.2  Wake  Flow  Region 

The  lower  wake  region  consists  of  two  subregions  separated  by  a 
line  of  zero  longitudinal  velocity  as  shown  in  Fig.  5.  The  upper  sub- 
region  consists  of  a  flow  characterized  by  the  dividing  streamline ‘ 
velocity  and  the  lower  subregion  represented  by  a  back  flow  velocity". 
A  linear  velocity  profile  is  assumed  for  the  upper  layer  while  a 
cosine  profile  is  assumed  for  the  back  flow.* 

The  mass  flow  in  the  upper  subregion  is 


fh  f1 

J  pudy  =  pdud6b(l  -  j- 


r) 


which  after  carrying  out  the  integration  for  isoenergetic  flows  gives, 

I 

>b 


pudy  = 


2JL 


d  d 


R<7  -  1)  Tq  2cd 


-£n(l  -  c2d) 


(2.12) 


For  the  ba/k  flow, 

/ 


/ 


p>  1 

J  pudy  -  pbubhb  J  £  ^  d{ 


or  / 


A 

/ 

J 

0 


pudv  =  —  p..  u  h  fl  -  I 
y  '  it  rb  b  b  \  b/j 


cos  ^  d^ 


0 


X-Cb 


2  2  „ 
cos  £ 


The  above  expression  may  be  integrated  directly  to  yield 


♦Accounting  for  the  wall  boundary  layer  does  not  change  the  basic 
equations  for  the  same  assumed  profile. 
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From  Eqs.  (2.12)  and  (2.13),  the  mass  balance  for  the  wake  region 
imposes  a  condition  that 
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In  addition,  if  a  straight  line  trajectory  is  assumed  for  the  iividing 
streamline,  the  geometry  of  the  wake  gives  the  following  relationships: 
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The  difference  in  pressure  across  the  wake  flow  may  be  calculated 
by  using  vertical  momentum  considerations.  The  control  volume  used 
for  this  analysis  is  shown  in  Fig.  6.  For  simplification,  it  is 
assumed  that  the  forward  flow  has  the  constant  pressure  p^  of  the 
dividing  streamline  while  the  back  flow  has  the  constant  wall  pressure 
Pw*  It  is  expected  that  the  normal  momentum  fluxes  associated  wj th  the 
transverse  velocity  component  are  small  and  tend  to  cancel  each  other. 
Thus,  they  are  neglected  in  this  analysis.  One  now  obtains 
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Upon  integrating  and  rearranging,  Eq.  (2.17)  becomes 
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We  now  consider  the  momentum  relationship  associated  with  the 
wake  flow.  The  momentum  associated  with  the  upper  subregion  is 
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which  after  integration  gives 
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For  the  lower  subregion,  the  momentum  associated  with  the  cosine 
profile  is 
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After  some  algebraic  manipulation,  Eq.  (2.20)  becomes 
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After  additional  rearranging  and  substitution  of  relationship  for 

dx  /'dx,  we  obtain 
w  ’ 


Using  Eqs.  (2.19)  and  (2.21)  and  normalizing,  the  above  equation 
becomes 


The  shear  stress  in  Eqs.  (2.7)  and  (2.;.. 5)  is  evaluated  from 
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which  can  be  rearranged  into 
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The  initial  conditions  for  this  recompressive  flow  are  provided 
for  frou  the  upstream  constant  pressure  mixing  region.  In  order  to 
achieve  a  smooth  joining  between  the  mixing  region  and  the  recomp^es- 
sion  region,  a  constant  pressure  laminar  mixing  analysis  was  performed 
by  assuming  a  velocity  profile  which  was  compatible  with  that  assumed 
for  the  reattachment  study. 

2.2.3  Laminar  Mixing  Analysis 

The  laminar  mixing  analysis  used  is  compatible  with  the  reattach- 
men*  analysis  but  applies  to  the  case  of  isoenergetic  flow  only.  The 
following  formulation  is  applicable  for  the  flow  with  a  finite  initial 
boundary  layer  thickress. 

Chapman  [9]  has  shown  from  a  rigorous  analysis  that  the  nondi- 
mensional  mixing  ordinate  for  similar  laminar  flow  can  be  represented 
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From  the  differential  equation  for  the  mixing  profile, 
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Chapman  obtained  the  solutions 


<pd  =  f '  CO)  =  0.5873 

and 
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These  equations  and  solutions  were  developed  by  assuming  that  the 
viscosity  was  a  function  of  temperature  only  and  additionally  that  no 
inicial  boundary  layer  thickness  was  present.  This  yields  a  single¬ 
valued  solution  for  the  dividing  streamline  velocity  of  0.5873. 

Further  research  has  shown  that  this  solution  is  valid,  however,  only 
for  large  distance  (x  »  1)  downstream  of  the  separation  point,  which 
is  not  the  case  for  realistic  base  flow  problems.  Additionally,  one 
recognizes  that  there  must  be  some  rotational  effects  when  the  flow 
field  interacts  with  the  upper  shear  layer. 

One  could  consider  the  rotational  flow  field  by  using  the  method 
of  characteristics  but  this  complicates  the  analysis  and  destroys  the 
simplicity  of  a  short  engineering  approach.  Dennison  and  Baum  [21] 
postulated  a  method  for  determining  the  dividing  streamline  velocity 
with  an  initial  boundary  layer  thickness.  They  uncoupled  the  momentum 
equation  from  the  energy  equation  for  a  finite  initial  profile  in  a 
laminar  free  shear  layer.  They  solved  the  resulting  equation  for 
velocity  and  shear  profiles  using  an  implicit  finite  difference  tech¬ 
nique.  The  initial  profile  as  the  flow  separates, corresponds  to  the 
Blasius  profile  in  transformed  coordinates.  Initially,  a  Howarth 
transformation  was  used  to  reduce  the  equations  to  the  incompressible 
form  and  for  convenience  of  calculation  another  transformation  to  the 
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required  for  the  reattachment  analysis  was  determined  from  a  logarithmic 
curve  fit  of  the  data  presented  in  this  figure. 


Equation  (2.26)  can  be  rewritten  for  the  upper  shear  ^.ayer  as 
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which  when  integrate 1  can  be  written  as 
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In  addition,  the  mixing  ordinate;  can  be  defined  as 


r,  = 


(2.35) 


which  is  the  same  as  that  for  turbulent  mixing  with  the  change  tp  the 

i  ■ 

appropriate  spread  rate • parameter  a.  The  equation  used  for  the  spread 
rate  parameter  is  that  used  by  Page  and  Dixon  [23],  i.e.. 


cr,  = 


(2.36) 


where  a  is  a  'eference  perturbation  velocity  factor,  u^/u^,  and  for 
this  analysis  is  taken  to  be  equal  to  2  in  accordance  with  the  analy¬ 
sis  of  Nash  [24] .  ' 
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From  the  exact  solution  for  laminar  mixing, 
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The  initial  values  required  for  the  problem  solution  are  obtained 
by  calculating  the  value  of  the  dividing  streamline  velocity  from  the 
curve  fit,  making  an  initial  guess  at  the  value  of  the  parameter 


5 a  '/Rex  /*  j  and  then  iterating  the  system  of  Eqs.  (2.33),  (2.34), 


and  (2.38)  until  the  system  converges  to  within  satisfactory  accuracy. 

From  this  analysis,  the  initial  valjes  of  &  ,  &,  ,  m, ,  and  the  mixing 

aba 

integrals  are  determined  for  use  in  the  reattachment  part  of  the 
analysis. 


2.3  Analysis  for  Redevelopment  Region 
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re^ttachment  point  is  carried  out  using  basically  the  same  techniques 
i  employed  for  the  reattachirent  analysis.  It  is  recognized  that  the 
characteristics  of  the  flow  will  be  the  realignment  of  the  external 
inviscid  flow  into  the  original  horizontal  direction.  The  initial 


conditions  required  can  be  estimated  from  previous  information  obtained 
at  the  reattachment  point.  Figure  8  is  a  representation  of  the 
geometry  and  the  modifications  required  to  establish  the  starting 
conditions. 

2.3.1  The  Flow  Condition  at  the  Initial  Section  of  Redevelopment 

Immediately  after  the  reattachment ,  the  viscous  flow  undergoes 
a  change  until  it  reaches  a  section  r‘  while  the  same  free-stream  con¬ 
dition  prevails.  A  continuous  rise  in  wall  static  pressure  has  been 
observed  and  it  can  only  be  made  possible  from  a  continuous  transfer 
of  mechanical  energy  to  the  dividing  streamline  (wall  streamline). 

Thus,  the  cytrvature  of  the  velocity  profile  at  the  wall  plays  a  pre¬ 
dominant  role  in  this  redevelopment  process.  The  assumed  velocity 
profile  for  the  reattaching  flow  does  possess  such  a  built-in  feature, 
namely  that  the  curvature  of  the  velocity  profile  at  the  wall  is  the 
largest  when  the  transfer  of  energy  is  most  needed  (otp/St;  i  d  =  0).  One 
may  thus  adopt  the  same  velocity  profile  for  the  entire  developing 
region  ((p^  =  0) , 

It  is  then  postulated  that  the  velocity  profile  slope  at  the 
new  section  r'  is  determined  from  the  condition  that  the  same  mass 
rate  of  flow  passes  through  both  sections  r  and  r*  as  shewn  in  Fig.  8. 

From  continuity  principles,  one  may  writs 
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The  geometric  configuration  gives  the  relationship 
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Thus,  Eq.  (2.39)  can  be  written 
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Using  the  assumed  velocity  distribution 
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at  the  section  r’,  one  may  iterate  on  the  slope  dqp/df ! 0  until  the 
relationship  given  in  Eq.  (2.40)  is  satisfied. 

The  normal  momentum  conservation  principle  will  be  subsequently 
applied  to  determine  the  pressure  level  at  r'.  One  recognizes  that 
the  momentum  flux  normal  to  the  wall  at  the  section  is  given  by 
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After  some  algebraic  manipulations  and  integration  of  the  last  term, 
the  previous  equation  becomes 
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(2.44) 


This  equation  allows  the  evaluation  of  the  wall  pressure,  p^/pg,  at 


section  rs . 

2.3.2.  Downstream  Redevelopment  Flow 

Once  the  flow  properties  at  section  j>’  have  been  established,  it 
is  simple  to  see  that  the  downstream  flow  redevelopment  involves  the 
continued  adjustment  of  the  free  stream  toward  the  horizontal  flow 
direction  and  the  analysis  suggested  for  the  upper  viscous  layer  in 
the  preceding  reattaching  flows  may  readily  be  adopted  for  this 
rehabilitation  process.  The  two  equations  already  derived  in  Section 
(2.2.1)  are  rewritten  for  the  new  system  of  coordinates  pertaining  to 
this  redeveloping  flow  region  as 
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where  =  u(c^)  -  u(ce)  and  is  negative  throughout  this  region. 

One  important  feature  which  is  associated  with  this  flow  redevel¬ 
opment  can  be  discussed  in  detail  as  follows.  During  thj.s  flow 
redevelopment,  the  free  stream  adjusts  itself  continuously  toward  the 
horizontal  flow  direction  thereby  inducing  an  appreciable  difference 
in  the  static  pressure  across  the  viscous  layer.  On  the  other  hand, 
the  low  energy  viscous  layer  close  to  the  wall  has  to  receive  an  ade¬ 
quate  amount  of  mechanical  energy  to  cope  with  the  rise  in  pressure 
and  this  transfer  is  properly  represented  by  the  curvature  of  the 
velocity  profile  at  the  wall.  It  would  seem  to  be  advantageous  to 
suggest  that  this  curvature  shall  be  proportional  to  the  normal  pres¬ 


sure  difference,  i.e., 
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again  with  the  value  of  k  evaluated  at  the  section  r'.  This  arrange¬ 
ment  would  insure  that  the  condition  of  a  completely  rehabilitated  or 


relaxed  state  is  reached  far  downstream. 


According  to  Eq.  (2.11),  the  normal  pressure  gradient  for  this 


region  can  be  rewritten  as 
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has  been  employed. 
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3 .  CALCULATION  METHOD  AND  RESULTS 

3.1  Calculation  Procedures 

For  a  specific  flow  condition  where  the  free-stream  Mach  number 

M  and  the  boundary  layer  approaching  the  base  are  given  (the  boundary 

00 

layer  can  usually  be  calculated  from  Re  ),  one  may  select  the  base 

i/ 

pressure  ratio  p^/p^  and  the  length  along  the  wake  boundary  where 
recompression  begins.  The  initial  values  required  to  calculate  the 
re  compression  process  may  be  determined  by  using  the  method  outlined 
in  Section  2.2.3.  One  now  starts  to  integrate  the  system  of  differen¬ 
tial  equations  for  the  reattachment  region.  At  each  increment  of 
length  along  the  dividing  streamline,  values  for  the  free-stream 
Crocco  number  ce  and  the  dividing  streamline  velocity  are  selected 
and  iterated  upon  until  the  system  of  equations  is  satisfied.  For 
each  set  of  these  values,  6a  may  be  found  from  Eq.  (2.2),  P^/Pg  from 
Eq.  (2.11),  from  Eq.  (2.15),  Pw/Pd  from  Eq.  (2.18),  hfe  from  Eq. 
(2.16),  and  c^  from  Eq.  (2.14).  Using  this  information,  one  obtains 
two  residues  from  Eqs.  (2.7)  and  (2.23).  The  values  for  cg  and  <pd  at 
this  location  are  iterated  upon  until  the  residues  become  negligible. 

The  above  calculations  are  continued  until  the  stagnation  point  is 
reached.  At  this  point  <pd  is  set  to  zero  and  the  wall  pressure  now 
becomes  the  stagnation  pressure  of  both  the  dividing  streamline  and  the 
back  flow.  Using  the  normal  momentum  relationship  given  by  Eq.  (2.18), 
the  correct  value  of  the  free-stream  Crocco  number  at  the  stagnation 


point  can  be  determined  when  the  wall  pressure  is  obtained  frc  n  the 
averaging  of  previously  determined  curves  for  pQj  and  pg^.  Once  again, 
two  residues  are  obtained  at  the  rear  stagnation  point.  These  are 
the  residues  for  the  whole  set  of  calculations.  One  now  iterates  upon 
the  initially  selected  values  for  base  pressure  ratio  P^/p^  and  the 
location  &  for  the  starting  of  the  recompression  process  until  the 
residues  at  the  stagnation  point  becomes  vanishingly  small.  This 
establishes  the  correct  flow  field  up  to  this  point. 

From  the  preceding  calculations  one  has  available  the  initial 
values  needed  for  calculating  the  flow  field  redevelopment  downstream 
of  the  reattachment  point.  Using  Eqs.  (2.44),  the  value  for  the  wall 
pressure  ratio  Pw/p^  is  determined  for  the  section  r'.  Then  Eqs. 
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In  addition,  the  lateral  shear  stress  terms  have  been  neglected  in  the 
estimation  of  the  pressure  differences  across  the  layers  since  their 
contribution  was  found  to  be  indeed  small. 

3.2  Results  of  Calculations 

Figure  9  shows  a  comparison  between  the  calculated  and  experi¬ 
mental  wall  pressure  ratio  Pw/P-^  as  a  function  of  downstream  distance 
from  the  step  for  a  free-stream  Mach  number  of  2.0  and  a  Reynolds 
number,  based  upon  the  characteristic  length  £  from  the  leading  edge 
of  the  plate  to  the  point  of  separation,  of  170,000  [25],  The  calcu¬ 
lated  base  pressure  value  of  0.680  compares  very  favorably  with  the 
experimental  value  of  0.670;  in  fact,  this  represents  an  error  of  only 
approximately  2  percent.  The  pressure  rise  in  the  recompression 
region  is  steeper  for  the  theoretical  solution  than  that  for  the 
experimental  data  but  the  overall  comparison  looks  satisfactory.  The 
corresponding  velocity  and  stagnation  pressure  ratios  for  this  case 
are  shown  in  Figs.  10  and  11. 

It  is  pertinent  to  mention  that  other  research  activities  on 
separated  flow  conducted  at  the  University  of  Illinois  [26]  indicated 
that  the  point  of  reattachment  clearly  exhibits  itself  as  a  saddle 
point  for  the  governing  system  of  equations.  This  fact  implies  that 
the  present  calculations  may  predict  accurately  the  base  pressure  but 
not  the  flow  properties  in  the  vicinity  of  the  point  of  reattachment. 
Indeed,  if  one  continued  the  calculations  for  the  redevelopment  region 
for  the  flow  case  presented  in  Fig.  9,  the  wall  pressure  ratio  would 
not  have  agreed  with  the  experimental  data  within  that  region.  Adjust¬ 
ments  in  the  thickness  of  the  viscous  layer  and  the  free-stream  Mach 
number  have  been  made  at  the  point  of  reattachment  such  that  the 


calculated  results  agreed  favorably  with  the  experimental  data.  This 
does  indicate  that  the  scheme  of  analysis  for  the  flow  rehabilitation 
downstream  of  the  reattachment  would  describe  the  flow  events  cor¬ 
rectly  if  the  conditions  at  the  point  of  reattachment  were  accurately 
determined . 

One  of  the  ways  to  assess  the  validity  of  a  particular  analytical 
method  is  to  assess  its  application  at  several  different  Mach  numbers. 
Some  difficulty  was  noted  when  the  proposed  method  was  applied  for 
flows  with  higher  Mach  numbers.  The  base  pressure  ratio  correlated 
very  well  with  experimental  data  [27]  over  a  wide  Reynolds  number 
range  at  a  Mach  number  of  2.5  as  shown  in  Fig.  12;  however,  the  wall 
pressure  ratio  rise  within  the  reattachmeut  region  did  not  agree  at 
all.  Again,  this  perhaps  shows  the  saddle  point  behavior  of  the 
reattachment  point.  Other  factors  which  could  have  caused  such  a  dis¬ 
crepancy  wili  be  discussed  later.  The  excellent  correlation,  however, 
for  the  base  pressure  ratio  over  the  large  Reynolds  number  range 
further  substantiates  the  fact  that  a  unique  solution  for  the  base 
pressure  is  obtained  from  the  solution  of  the  integral  equations. 

Figure  13  shows  the  calculated  wall  shear  stress  downstream  of 
the  reattachment  point.  The  distribution  appears  to  follow  the 
expected  trend  beginning  with  zero  at  the  stagnation  point,  rising  to 
a  plateau  level  thereafter.  After  the  flow  reaches  the  fully  rehabil- 
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back  step  is  disregarded,  is  also  shown  i.  Fig.  13.  This  curve  serves 
to  indicate  the  correct  level  of  the  plateau  obtained  from  the  shear 
stress  calculations. 


4.  DISCUSSION  AND  CONCLUSIONS 


The  reattachment  and  subsequent  redevelopment  for  laminar  super¬ 
sonic  flow  over  a  back  step  can  be  calculated  using  the  described 
procedure.  The  calculated  results  compare  favorably  for  one  particular 
Mach  number  while  yielding  a  generally  unsatisfactory  pressure  distri¬ 
bution  at  a  higher  Mach  number.  Several  reasons  for  this  disagreement 
are  apparent.  First  of  all,  it  is  difficult  to  evaluate  or  to  com¬ 
pletely  satisfy  one's  self  that  the  experimental  data  are  for  com¬ 
pletely  laminar  flow  even  though  the  Reynolds  numbers  are  low.  In 
fact,  if  the  flow  is  transitional,  one  would  expect  the  pressure 
gradient  and  the  reattachment  point  to  occur  closer  to  the  base  than 

that  for  laminar  flow  and  this  is  the  case  for  the  data  used  at  Mach 

* 

number  2.5.  Additionally,  one  may  expect  that  as  the  Mach  number  is 
increased,  disregarding  the  rotational  flow  field  interacting  with  the 
upper  viscous  layer  might  have  caused  a  sizable  error  In  the  estimation 
of  the  initial  conditions  required  for  the  reattachment  calculations. 
This  effect  was  recognized  in  the  latter  part  of  this  investigation. 

It  is  suggested  that  future  research  studies  should  take  this  effect 
into  consideration. 

The  method  chosen  for  the  redevelopment  region  has  been  shown  to 
yield  satisfactory  results.  Accounting  for  the  wall  shear  stress  made 
little  difference  in  the  wall  pressure  ratios.  This  probably  was  due 
to  the  exceedingly  small  magnitude  of  the  calculated  shear  stress 


values. 
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Additionally,  it*was  also  determined  through  sample  calculations 
that  the  inclusion  of  the  lateral  shear  stress  in  the  estimation  of 
the  pressure  difference  across  the  viscous  layer  made  essentially  no 

i 

change  in  the  magnitude  of  the  bass  pressure  and  the  distribution  of 
the  wall  pressure  ratio.  This  fully  justifies  the  original  approach 
to  disregard  their  contribution.  . 

Finally,  it  should  be  .noted  that  this  suggested  flow  model  can  be 
used  for  other  separated  flow  problems.  Similar  models  can  be  devel¬ 
oped  for  low  speed  flow  and  the  analysis  can  also  be  extended  -to 
axisymmetrical  configurations.  :  > 
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Figure  9 


Figure  10. 


.  Theoretical  and  Experimental  Wall  Pressure  Distribution 
for  M  =  2.0,  R  =  170,000 

& 


Figure  13.  Variation  of  Wall  Shear  Stress  in  Redevelopment  Region 

for  M.  =  2.0,  R  =  170,000 
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